FOURIER SERIES AND APPROXIMATION ON HEXAGONAL 
AND TRIANGULAR DOMAINS 



YUAN XU 



Abstract. Several problems on Fourier series and trigonometric approxima- 
tion on a hexagon and a triangle are studied. The results include Abel and 
Cesaro summability of Fourier series, degree of approximation and best approx- 
imation by trigonometric functions, both direct and inverse theorems. One of 
the objective of this study is to demonstrate that Fourier series on spectral sets 
enjoy a rich structure that allow an extensive theory for Fourier expansions 
and approximation. 



1. Introduction 

A theorem of Fuglede [6] states that a set tiles R" by lattice translation if and 
only if it has an orthonormal basis of exponentials e*^"^^^ with a in the dual lattice. 
Such a set is called a spectral set. The theorem suggests that one can study Fourier 
series and approximation on a spectral set. For the simplest spectral sets, cubes 
in M'', we are in the familiar territory of classical (multiple) Fourier series. For 
other spectral sets, such a study has mostly remained at the structural property 
in the level and it seems to have attracted little attention among researchers in 
approximation theory. 

Besides the usual rectangular domain, the simplest spectral set is a regular 
hexagon on the plane, which has been studied in connection with Fourier anal- 
ysis in [H [13] . Recently in [10] , discrete Fourier analysis on lattices was developed 
and the case of hexagon lattice was studied in detail; in particular, Lagrange in- 
terpolation and cubature formulas by trigonometric functions on a regular hexagon 
and on an equilateral triangle were studied. Here we follow the set up in [TU] to 
study the summability of Fourier series and approximation. The purpose of this 
paper is to show, using the hexagonal domain as an example, that Fourier series 
on a spectral set has a rich structure that permits an extensive theory of Fourier 
expansions and approximation. It is our hope that this work may stimulate further 
studies in this area. 

It should be mentioned that, in response to a problem on construction and 
analysis of hexagonal optical elements, orthogonal polynomials on the hexagon were 
studied in [5] in which a method of generating an orthogonal polynomial basis was 
developed. We study orthogonal expansion and approximation by trigonometric 
functions on the hexagon domain. 
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In comparison to the usual Fourier series for periodic functions in both variables 
on the plane, the periodicity of the Fourier series on a hexagonal domain is defined 
in terms of the hexagon lattice, which has the symmetry of the reflection group A2 
with reflections along the edges of the hexagon. The functions that we consider are 
periodic under the translation of hexagonal lattice. It turns out ([10l[T3]) that it is 
convenient to use homogenous coordinates that satisfy ti + i2 + ^3 = in M.^ rather 
than coordinates in . Using homogenous coordinates allows us to treat the three 
directions equally and it reveals symmetry in various formulas that are not obvious 
in coordinates. As we shall show below, many results and formulas resemble 
closely to those of Fourier analysis and approximation on 27r periodic functions. 

Fourier analysis on the hexagonal domain can be approached from several direc- 
tions. Orthogonal exponentials on the hexagonal domain are related to trigonomet- 
ric functions on an equilateral triangle, upon considering symmetric exponentials 
on the hexagonal domain. These trigonometric functions arise from solutions of 
Laplacian on the equilateral triangle, as seen in |^ and developed extensively in 
(13. ,14 , and they are closely related to orthogonal algebraic polynomials on the do- 
main bounded by Steiner's hyercycloid [SlITO], much as Chebyshev polynomials arise 
from exponentials. In fact, the trigonometric functions arise from the exponentials 
by symmetry are called generalized cosine functions in [101 113j , and there are also 
generalized sine functions that are anti-symmetric. Our results on the hexagonal 
domain can be easily translated to results in terms of generalized cosines. Our 
results on summability can also be translated to orthogonal expansions of algebraic 
polynomials on the domain bounded by hypercycloid, but the same cannot be said 
on our results on best approximation. In fact, just like the case of best approxi- 
mation by polynomial on the interval, the approximation should be better at the 
boundary for polynomial approximation on the hypercycloid domain. For exam- 
ple, our Bernstein type inequality (Theorem 4.8) can be translated into a Markov 
type inequality for algebraic polynomials. A modulus of smoothness will need to 
be defined to take into account of the boundary effect of the hypercycloid domain, 
which is not trivial and will not be considered in this paper. Some of our results, 
especially those on the best approximation, can be extended to higher dimensions. 
We choose to stay on the hexagonal domain to keep an uniformity of the paper and 
to stay away from overwhelming notations. 

The paper is organized as follows. Definitions and background materials will be 
given in Section 2. In Section 3 we study the Abel summability, aka Poisson integral, 
and Cesaro (C, S) means of the Fourer series on the hexagon, where several compact 
formulas for the kernel functions will be deduced. One interesting result shows that 
the (C, 2) means are nonnegative, akin to the Fejer means for the classical Fourier 
series. In Section 4 we study best approximation by trigonometric functions on 
the hexagonal domain and establish both direct and inverse theorems in terms of 
a modulus of smoothness. 



Below we briefly sum up what we need on Fourier analysis on hexagonal domain. 
We refer to [10] for further details. The hexagonal lattice is given by HZ'^ , where 
the matrix H and the spectral set fin are given by 



2. Fourier series on the regular hexagon 
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respectively. The reason that ^Ih contains only half of its boundary is given in ^10]. 
We will use homogeneous coordinates (<i,i2,i3) that satisfies ti + t2 + ts = for 
which the hexagonal domain Qh becomes 

n = {iti,t2.,t3) eR^ : -1 <ii,i2,-i3 < 1; h + t2 + t3^0} , 

which is the intersection of the plane ti + ^3 = with the cube [—1, 1]"^, as seen 
in Figure 1. 



Figure 1. Regular hexagon in and in R^. 
The relation between {xi, X2) G and t G f2 is given by 

(2-1) ^^~Y^~' t2=X2, *3 = -Y —■ 

For convenience, we adopt the convention of using bold letters, such as t, to denote 
points in the space 

R% ~ {t = {ti,t2,h) e : ti + i2 + ^3 - 0}. 

In other words, bold letters such as t stand for homogeneous coordinates. If we 
treat a; g and t e as column vectors, then it follows from (|2.ip that 

(2.2) X - \H{h -hM- hf = \H{2h +t2,h + 2t2f' 

upon using the fact that ti + 12 + ^3 = 0. Computing the Jacobian of the change of 
variables shows that dx = ^^dtidt2- 

A function / is called periodic with respect to the hexagonal lattice, if 

f{x)^ f{x + Hk), k£l?. 

We call such a function _ff-periodic. In homogeneous coordinates, x = y (mod H) 
becomes, as easily seen using ()2.2p . t = s mod 3, where we define 

t = s mod 3 ti ^ si = t2 — S2 = is — S3 mod 3. 

Thus, a function /(t) is iJ-periodic if /(t) = /(t + j) whenever j = mod 3. If / 
is i7-pcriodic, then it can be verified directly that 

(2.3) / /(t + s)dt== / /(t)dt, seR^H- 
Jn Jn 
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We define the inner product on the hexagonal domain by 

(/,.9)h := TTT-T / f{xi,X2)g{xi,X2)dxidx2 ^ -r^ / /(t)5r(t)c?t, 
I^'hI JnH l^'l 

where \n\ denote the area of fl. Furthermore, let Z,^ := 1? n R|^. We define 

0j(t) :=e^J-*, i^I?H, t^W\j. 

These exponential functions are orthogonal with respect to (/, .g)^. 

Theorem 2.1. 6,. for k,j e Z|^, (0k, 0j) = (^kj- Moreover, the set : j G Zf^} 
zs an orthonormal basis of L^{^1). 

It is easy to see that 0j are iJ-periodic functions. We consider a special collection 
of them indexed by j inside a hexagon. Define 



■.-n< ji, j2, < n} and J„ ~ H„ \ H„_ 



1- 



Notice that j e ]HI„ satisfies ji + j2 + = 0, so that Mn contains all integer points 
inside the hexagon n fi, whereas J[„ contains exactly those integer points in that 
are on the boundary of nft, or n-th hexagonal line. We then define 

(2.4) H„ :-span{0j :j eH„}. 

It follows that dim Tin — + 3n + 1. As we shall see below, the class Tin shares 
many properties of the class of trigonometric polynomials of one variable. As a 
result, we shall call functions in Tin trigonometric polynomials over Q. We will 
study the best approximation by trigonometric polynomials in 7i„ in Section 4. 

By Theorem 2.1, the standard Hilbert space theory shows that an i?-periodic 
function / G L^(ri) can be expanded into a Fourier series 

(2.5) f^J^y^ inL^{Q), where fy.= ±. f f{s)e-'^ ds. 

We consider the n-tli hexagonal Fourier partial sum defined by 

(2.6) Snfix) := J2 fM^) = ^ / /(t - s)i?„(s)ds, 

where the second equal sign follows from (j2.3p with the kernel £)„ defined by 

jeH„ 

The kernel is an analogue of the Dirichlet kernel for the ordinary Fourier series. It 
enjoys a compact formula given by [13] (see also [TU]) that 

(2.7) i?„(t) = e„(t)-e„_i(t), 

where 

(2.8) e„(t) = 3 3 3 



sm ^ 3 ^ sm ^ sm 



This last formula is our starting point for studying summability of Fourier series in 
the following section. 



"'^Here our notation differs from tliose in .lOi. Our EI,i and 7i„ are in fact H* and Ti.* tliere. 
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We will also need to use the Poisson summation formula associated with the 
hexagonal lattice. This formula takes the form (see, for example, [1] \7\ ITO]) 

(2.9) J2 /(^ + = JTTm E 7(^-'^A)e^-'^"^"^^ 

where for / G L^(]R^) the Fourier transform / and its inverse are defined by 

m = / /(x)e-2-«-da; and f{x) = / /(^e^-^-dx, 

and the formula holds under the usual assumption on the convergence of the series 
in both sides. Define Fourier transform in homogeneous coordinates by /(t) := 
f{H^'^'t) with t = (ti,i2) as a column vector. Using (|2.2p . it is easy to see that / 
and its inverse become, in homogeneous coordinates, 

(2.10) 7(s) = / /(t)e-^^Vt and /(t) = / 7(s)e^^*ds. 

Next we reformulate Poisson summation formula in homogenous coordinates. The 
left hand side of (|2.9p is an //-periodic function over hexagonal lattice, which 
becomes the summation of /(t + 3j) over all j e Z|j as x = y mod H be- 
comes t = s mod 3. For the right hand side, using the fact that {ki, k2)H^^x = 
i(fci, /C2)(ii -hM- hT = • t by (I22D, we obtain 



~2irfc"//-^x , _ 2V3 [ „,,^ _27rkt 2\/3- 



f{H-'^k) = / /(x)e-2-^- " -dx=^ /(t)e-2-k t^t ^ ^/(k). 



Consequently, we conclude that the Poisson summation formula in (j2.9p becomes, 
in homogeneous coordinates, 

(2.11) E /(t + 3k) = i ^ 7(j)e^J *. 

Throughout of this paper we will reserve the letter c for a generic constant, 
whose value may change from line to line. By A ^ B we mean that there are two 
constants c and c' such that cA < B < c'A. 

3. SUMMABILITY OF FOURIER SERIES ON HEXAGON 



We consider hexagonal summability of Fourier series (|2.5P : that is, we write the 
Fourier series p.Sp as blocks whose indices are grouped according to J„: 

00 

/(t)-EE7'/>j(t). 

From now on we call such a series hexagonal Fourier series. Its n-th partial sum is 
exactly 5„/ given in (|2.6p . 

3.1. Abel summability. We consider Poisson integral Prf{t) of the hexagon Fourier 
series, which is defined by 

r^r:if "Jim 
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where P(r;t) denotes the Poisson kernel 

P(r;t) :=^^0j(t)r", < r < 1. 

n=OjGJI„ 

Just as in the classical Fourier series, the kernel is nonnegative and enjoys a closed 
expression. 

Proposition 3.1. (1) The Poisson kernel P{r;t) is nonnegative for a/Z t G and 

(2) Let q{r, i) = 1 - 2r cos i + . Then 

(1 - rf{l - r^) 



P(r;t) 
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r(l — r)^ r(l — r)^ 



r 



(1-r) 



2 



Proof. The fact that the integral of P{r; t) is 1 follows from the definition and the 
orthogonality of whereas that P{r; t) > is an immediate consequence of the 

compact formula in the part (2). 

To prove the part (2), we start from the compact formula of Z)„(t), from which 
follows readily that 

oo oo 

P(r; t) = (1 - r) ^ i?„(t)r" = (1 - rf ^ e„(t)r". 

n=0 n=0 

li ti + t2 + t^ — then it is easy to verify that 

sin 2ti + sin 2t2 + sin 2^3 = — 4 sin ti sin t2 sin ^3 , 

(3.1) 

cos 2fi + cos 2^2 + cos 2f3 = 4 cos ti cos t2 cos f3 — 1. 

Using the first equation in (|3.H) and the fact that 

00 

sin(n + l)s = ^ 

we conclude then 



2r cos s + r^ 

n—O 



1 + 2r + 2r'^ + r** - 2r'^ 


27r(ti-t2) , 27r(ti-t3) , 27r(t2-t3)" 
COS ^ + COS ^ + COS ^ 









27r(t3-ti) 
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Putting the three terms together and simpUfy the numerator, we conclude, after a 
tedious computation, that 

oo 
n=0 

The numerator of the right hand side can be written as 

(1 - r)(l -r^)+r [(r, HlillzM) + q + q 

from which the stated compact formula follows. □ 

Next we consider the convergence of Prf as r 1 — . For the classical Fourier 
series, if Prf converges to / as r ^ 1— then the series is called Abel summable. 

Theorem 3.2. // / is an H -periodic function, bounded on f2 and continuous at 
t e ^l° , then Prf{t) converges to /(t) as r i— > 1—. Furthermore, if f is continuous 
on fl, then Prf converges uniformly on Q as r I— . 

Proof. Since / is continuous at t e For any e > 0, choose 6 > such that 
|/(t - s) - /(t)| < e whenever s e :^ {s e Q. : \si\ < (5, 1 < i < 3}. 
Since P(r; s) has a unit integral, it follows that 

l^r./(t) - /(t)| < e [ P{r; s)ds + 2 / |/(t - s) - /(t)|P(r; s)ds 

<e + 2\\f\\^. f P{r;s)ds. 
jn\ns 

Thus, it suffices to show that the last integral goes to as r i-^ 1+. Since q{r,t) = 
(1 — r)^ + 2r sin^ | > (1 — r)^, the closed formula of the Poisson kernel shows that 
Pr{t) is bounded by 

, , 2(1 -rf 2(1 -r)2 



, 2(1 -r)^ 

q [r, H£(£|^) q Mii^ 

Clearly we only need to consider one of the three terms in the right hand side, say 
the second one, so that the essential task is reduced to show that, as q{r,t) is even 
in t, 

(3.2) lim / — ^ — ^ds = 0. 



Setting P{r;t) := ^^yj^ for — tt < t < tt. Then P{r;t) is the Poisson kernel for the 
classical Fourier series. It is known [16] that 



(3.3) / P{r;t)dt = l and 0<P(r;t)<c 

J —TV 



1 -r 
1^' 
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where c is a constant independent of r and t. Evidently, the kernel in the left hand 
side of (13.21) can be written as 

^ ^ ^ p f „. 2£(ti-i3) \ p / 27r(f2-t3) \ 

Recall that + t2 + ^3 = for t = (ii, i2, ^3) G It is easy to see that t e 
means that ~1 < ti,t2,ti — t2 < I- Let 

, . ^1—^3 2ti+t2 t2 — ti+2t2 

(3.4) ,,:=^ = ^_, ,,:=__ = ^_. 

A simple geometric consideration shows that the image of the domain under 
the affine mapping (ii,t2) i~> (si,S2) in (|3.4p contains the square [—(5/6,(5/6]^. 
Consequently, the image of \ ll^ is a subset of [—1,1]^ \ [—5/6,6/6]'^. Hence, 
changing variables (ti,t2) ^ (si,S2), we obtain by (|3.3p that 



3 

(1 + r)2 y[„l,l]2\[_5/g ,5/g] 



dt 



^ ,1 , „^2 / ^ ('^' 27rsi) P (r, 27rs2)ds 



< 3 / P{r,2T:s)ds < c 

Jl>lsl>5/6 



1 - r 



'l>|s|>(5/6 ^ 

which converges to zero when r ^ 1—. This proves (|3.2p and the convergence 
of Prf[t) to /(t). Clearly it also proves the uniform convergence of P^/ to / S 

cip.). □ 

3.2. Cesaro summability. Let us denote by Sf^f and the Cesaro (C, (5) means 
of the hexagonal Fourier series and its kernel, respectively. Then 

^n/(t) := / /(s)x(^)(t6-s)ds, 
Jn 

where 

" fc=o jej„ ^ 

It is evident that the case (5 = corresponds to Snf and _D„(t), respectively. By 
and (H^, the (C, 1) kernel is given by 



if«(t) = le„(t)^i 



1 sin ("+^»*^-*-)- sin sin 



sm ^ 3 ' sm ^ sm ^ 3 " 



For the classical Fourier series in one variable, the (C, 1) kernel is the well-known 
Fejer kernel, which is given by 

and is nonnegative in particular. For the hexagonal Fourier series, it turns out that 
(2) 

Kn is nonnegative. 
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Lemma 3.3. For n > 0, 



A„(t) := cos ^ sin + cos ^ sin ^"+^^^^-^-*^^" + cos ^ sin 

B„(t) sin sin (»+2)(«2-t3)^ ^ ^ ^^^^ gin ("+2)fo-t2)^ ^ 

Proof. Using p.ip and the elementary formula 



fc=0 



sm[2(fc + 1).] = BM(n + l).]sin[(n + 2).] 



we obtain 

-4^„(t) 



:e,(t) = 



fc=o I sm ^ ^ ' 1 I sm 1 I sm ^ ^ ^ 1 

where 

E^{t) :^sin(^i±iKii^sinMI^i^ 



sin 


(t2 


-t3)7I- 

3 


sin 


(ts- 


-ti)7r 
3 


sin 


(tl 


-t2)7r 
3 


sin 


(t3- 


-ti)ir 
3 


sin 


(tl 


-t2)7r 
3 


sin 


(t2- 


-t3)7r 
3 



3 3 

+ Sin sin 

+ sin sin ("+^»*^-*^)- 

From here, the difficulty lies in identifying that the numerator is a sum of squares. 
Once the form is recognized, the verification that — 4_E„(t) ~ j^{{An{tY + -B„(t)^) 
is a straightforward, though tedious, exercise. □ 

An immediate consequence of the above lemma is the following: 

Theorem 3.4. The (C, 2) means of the Fourier expansion with respect to the 
hexagon domain is a positive linear opreator. 

As a comparison, let us mention that for the usual Fourier series on the tours, if 
the partial sums are defined with respect to £°° ball, that is, if the Dirichlet kernel 
is defined as 

then it is proved in [2] that the corresponding (C, 3) means are nonnegative and 
(5 = 3 is sharp. In fact, the results in 2] are established for the partial sums defined 
with respect to the £^ ball for d-dimensional torus. For d = 2, it is easy to see that 
the result on £^ ball implies the result on ball. 

If the (C, S) means converge to /, then so is (C, S') means for S' > S. The posi- 
tivity of the kernel shows immediately that (C, 2) means of the hexagonal Fourier 
series converges. It turns out that the (C, 1) means is enough. 

Theorem 3.5. /// G C{Q), then the (C, 1) means Sn^ f converge uniformly to f 
in VL. 
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Proof. A standard argument shows that it suffices to prove that Sn"^ is a bounded 
operator, which amounts to show that 

|e„(t)| dt < cn 

Since 0n(t) is a symmetric function in ^1,^2,^3, we only need to consider the integral 
over the triangle 

A := {t e Mh : < ii, i2, -ts < 1} = {(^1,^2) : <i > 0, i2 > 0, ti+t2< f }, 

which is one of the six triangle in fl (see Figure 1). Let si, S2 be defined as in (|3.4p 
and let A* denote the image of A in (si,S2) plane. Then 

A = {(si, S2) : < si < 2s2, < S2 < 2si, si + S2 < 1} 
and it follows, as the Jacobian of the change of variables is equal to 3, that 
sin((7i + l)7rsi) sin((n + l)7rS2) sin((7i + l)7r(si — S2)) 



dsids2 



sin(7rsi) sin(7rs2) sin(7r(si + S2)) 
Since the integrant in the right hand side is clearly a symmetric function of (si, S2), 
it is equal to twice of the integral over half of the A, say over 

A* = {(si, S2) e A : si < S2} = {(si, S2) : si < S2 < 2si, si + S2 < !}• 
Making another change of variables si — {ui — U2)/2 and S2 — {ui + U2)/2, the 
domain A* becomes F := {(ui,M2) : < M2 < ui/3,0 < mi < 1} and we conclude 
that 



duidu2 i j \^*n{u)\du 



iui-\-U2)'^ • (ui—U2)7r • U9TT 

sm ^ ^ ' sm 2 si'^ 2 
To estimate the last integral, we partition F as F Fi U F2 U F3 and consider the 
three cases separately. 

Case 1. Fi = {u G F : ui < 3/n}. Using the fact that | smnt/ sint| < n, we obtain 

r 3j^2 
\Ql{u)\du < {n + if / duidu2^ [n+lf— <cn. 

Case 2. F2 = {u G F : 3/n < mi, U2 < In this region we have U1 — U2 > 2ui/3 

and ui + U2 > ui. Hence, upon using sinu > {2/t:)u, we obtain 



3/ri 



sm 2 2 



1 1 

dui < I —dui < cn. 

3/n "1 



Consequently, it follows that 

"l/n /.I 

\e*ju)\du 





sin 


Jo 


sinifl^ 



< cn 



\e*„{u)\du < cn 

J 3/71 

upon using | sinnu/ sinu| < n again. 

Case 3. F3 = {tt G F : 3/n < u_ U2 > In this region we have ui > 3u2, which 

implies that U1 — U2 > (2/3)ui and U1+U2 >ui. Thus, using sin m > (2/tt)u again, 
we obtain 



' 3m2 



gjjj ("+l)("l+"2)7r {n+l){ui-U2)7T 



sm sm -^^ 2 — ^ 



1 , c 
aui < c / —^dui < — . 

3U2 ^1 ^2 
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Consequently, using sinu2 > (2/7r)u2, we conclude that 

/■I/" /•! /-i/s 1 

\e*Ju)\du^ / \ei{u)\du<c —du2<cn. 

Putting these estimates together completes the proof. □ 

This theorem shows that (C, S) summability of the hexagonal Fourier series be- 
haviors just like that of classical Fourier series. In particular, we naturally conjec- 
ture that the (C, 6) means of the hexagonal Fourier series should converge if 6 > 0. 
This condition is sharp as (5 = corresponds to 5„/ whose norm is known to be in 
the order of (logn)^ ([HI [13]). 

4. Best Approximation on Hexagon 

For 1 < p < oo we define space to be the space of Lebesgue integrable 
i/-periodic functions on Q with the norm 

i/p 



jjfit)\Pdt^ ^ i<p<oo. 



and we assume that is C{n) when p = oo, with the uniform norm on fl. For 
/ e LP, we define the error of best approximation to / from Tin by 

EM)p :=^infj|/-5||,. 
We shall prove the direct and inverse theorems using a modulus of smoothness. 

4.1. A simple fact. We start with an observation that En{f) can be related to 
the error of best approximation by trigonometric polynomials of two variables on 
[—1,1]^. To see this, let us denote by Tn the space of trigonometric polynomials of 
two variables of degree n in each variable, whose elements are of the form 

n n 

T{u) =^Y^ l,^^2Mkiui+k2U2) _ 
ki=Q k2=0 

Furthermore, if / G C([— 1, 1]^) and it is 2tt periodic in both variables, then define 
^n(/):=inf max J/H-r(u)|. 

Proposition 4.1. Let f be H -periodic and continuous over the regular hexagon 
fin- Assume that f* is a continuous extension of f on [—1,1]^. Then 

<fLf ](/*)■ 

Proof. We again work with homogeneous coordinates. Using the fact that ti+t2 + 
is = and ji + j2 + js = 0, we can write 

0j(t) = e^j-* = e2^«0i«i+J2«2)^ ^ 2*1^*2 ^ ^ tj^_ 

Clearly t G fl implies that (si,S2) G fi* {s : —1 < si,S2,si + S2 < 1}. Fur- 
thermore, j e IHI„ implies that —n < ji,j2,ji +32 < n. Consequently, we have 
that 

Suit) cj'^j(t) = E c,„,„_,,_,,e2"to^i+^-^^^) r„(s). 

-n<ji-\-j2<n 
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Let g{s) = /(t) = /(2si — S2,2s2 — si, — si — S2) and let its continuous extension 
to [-1, 1]2 be 5*. It follows that 

\\f-Sn\\=m^\g{s)-Tn{s)\ < max - T„(s)|. 

sGf2* — l<si,S2<l 

Taking minimum over all Sn € translates to taking minimal of all cj, conse- 
quently we conclude that 



En{f) < min max 

CJ — 1<S1,S2<1 



< min max 

Cj -1<S1,S2<1 



E 



.p27ri(ji«i+j2S2) 



CiC 



-n<ji+j2<n 



E 

-LfJ<ji.j2<LfJ 



j,jg27ri(jisi+j2S2) 



If we work with / defined on Qh, then g* becomes /* in the statement of the 
proposition. □ 

For smooth functions, this result can be used to derive the convergence order of 
En{f). The procedure of the proof, however, is clearly only one direction. Below 
we prove both direct and inverse theorems using a modulus of smoothness. 

4.2. Modulus of smoothness. On the set O, we can define a modulus of smooth- 
ness by following precisely the definition for periodic functions on the real line. 
Thus, we define for r e No, 

At/(x) = /(x + 1) - /(x), A[/(x) = At A[- V(x). 

It is well known that 



At/(x)=^(-l)'=Q/(x + fct). 

fc=0 ^ ' 



For t e Rh, let ||t|| := {t{ the usual Euclidean norm of (ii,i2)- We can 

also take other norm of (ti,t2) instead, for example, ||t||oo := max{|ti|, \t2\\- Since 
ti+t2+h = 0, we have evidently ||t||oo < max{|ti|, Itg], lis]} < 2||t||oo- The 
modulus of smoothness of an i7-periodic function / is then defined as 

Wr(/;/i)p:= sup ||A;/||j„ 1<p<oo. 
\\t\\<h 

Proposition 4.2. The modulus of smoothness satisfies the following properties: 

(1) For A > 0, ujr{f; \h)p < (1 + A)'^w^(/; h)p. 

(2) Let di denote the partial derivative with respect to ti and d'' = 81^82^8^^ 
for k = {ki,k2, ks). Then 



ki+k2+h3=r 



kt\k2\kiV 



Proof. The proof of (1) follows exactly as in the proof of one variable. For part (2), 
it is easy to show by induction that 



At7(x) 



/ 5«i • • • du. 



/(x + U\t -|- . . . + Urt)dU\ . . .Ur- 
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The integrant of the right hand side is easily seen, by another induction, to be 

ki+k2+k3=r 

from which the stated result follows from (|2.3p and the fact that Ht*^!! < /i'*^! = 
h\ □ 

4.3. Direct theorem. For the proof of the direct theorem, we use an analogue of 
Jackson integral. Let r be a positive integer. We consider the kernel 

K„,r{t) := Xn.r [6„(t)]^'' , whcrc / Kn,ri't)dt = 1. 
Since n~^8„(t) is the (C, 1) kernel of the Fourier series, we see that 8„ £ Hn and, 

thus, Kn,r G Hrn- 

Lemma 4.3. For e N and v <2r — 2, 

I \\tW'Kn,r(t)dt < cn-''. 
Jn 

Proof. First we estimate the constant Xn.r- We claim that 

(4.1) (A„,,)-' - / [e„(t)]''^dt^n- 

Jn 

We derive the lower bound first. As in the proof of Thcorem l3.51 we change variables 
from t g to (si, S2) G ^* , then use symmetry and change variables to (ui, U2) G F. 
The result is that 

/ [e„(t)]2'-dt = 3 / [e:(u)]^'^du>3 / [e:(u)]^'-du, 

where we choose F* — {{ui, U2) : i^^^^i-^ < "2 < ui/'i < ■^^^:;^:ppj} , which is a subset 
of F and its area is in the order of For it S F*, we have sin(rt + — U2) > 
sin7r/16, sin(n + l)(ui + U2) > sin7r/8, and sin(n + l)u2 > sin7r/32; furthermore, 
sin(Mi — U2) < sin57r/(16n) < 57r/(16?T.), sin(wi + U2) < sin7r/(2n) < 7r/(2n), and 
sinu2 < sin7r/(8n) < n/{8n). Consequently, we conclude that 



6r+2 



[ [e„(t)]^''dt > cn^'' [ du^ 



This proves the lower bound of (|4.ip . The upper bound will follow as the special 
case V — Q oi the estimate of the integral /^'"^ below. 
We now estimate the integral 



:= / |tr[e„(t)]^'-dt. 



Again we follow the proof of Theorem p.Sp and make a change of variables from t e 
to (wi, U2) G F. The change of variables shows that ^i = ^ [2(ui — M2) — [ui + U2)] 
and t2 = ^ [2(mi + U2) ~ ("i ~ ""2)], which implies that 

||t||oo = max{|ii|, \t2\} < 5 max{|ui - U2I, |ui + U2I} • 
Consequently, we end up with 

< c [ max{|«i - U2\, \ui + «2|}' mu)\^'du. 
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Since max{|a|, |5|} < \a\ + we can replace the max{...} term in the integrant 
by the sum of the two terms. The fact that < 2r — 2 shows that we can cancel 
\ui — U2\'^, or |ui + with the denominator of 8* . After this cancellation, the 
integral can be estimated by considering three cases as in the proof of Theorem l3.5l 
In fact, the proof follows almost verbatim. For example, in the case 2, we end up, 
using ui — U3 > 2ui/3 and ui + U2 > ui, that 



3/n 



{n+l){ui+U2)-!T 
2 

• {Ui+U2)1T 



-I 2 r 



{n+l){ui~U2)-n 



< 



sm 

1 

3/?i 



{ui—U2)Tr 



2r-v 



4r — 1/ 



dui < cn 



dui 



4i — u — 1 



Consequently, it follows that 

|ui-u2ne;(u)|du = 



1/ti 



< cn 



3/n 



\ui-U2V\Ql{u)\du 
2r 





(n-\-l)u2TT 

sm 2 


'0 





du < cn'^"-"-'^. 



upon using | sin nit/ sin m | < n. The other two cases can be handled similarly. As 
a result, we conclude that I^'^ < cn^^^P^^. The case p = gives the lower bound 
estimate of (|4.ip . The desired estimate is over the quantity Xn,rln^ and follows 
from our estimates. □ 

Using the kernel if„,r we can now prove a Jackson estimate: 

Theorem 4.4. For I < p < 00 and for each r — 1,2, . . there is a constant Cr 
such that if f £ then 

En{f)p < CrUJrif, ^)p. n=l,2, .... 

Proof. As in the proof of classical Jackson estimate for trigonometric polynomials 
on [0, 27r], we consider the following operator 

Fr/(x) := / J„,p(t)^(-l)'=-i('';V(x + fct)dt, 

where J„,p(t) = i^„.,p(t) with n* = [^J + 1, and p > {r + 2)/2. Evidently, 
Jn,p{—'t) — Jn.pit)- Using the fact that J„_p e 7i„, we see that Fn^pf can be 
written as a linear combination of 



(4.2) 



/ /(x + M)(/)j(t)dt, ieHn, fc = i, 
Jn 



As / is iJ-periodic, so is /(x + A:t) as a function of t. Let F^ — J2jeH„, '^i^i denote 
the (C, 1) means of the Fourier series of / over fi. Then Fm converges to / uniformly 
on n. If j ^ —kl for some 1 G Hn then, using the fact that <^i(x + A;t) = 0i(x)</)fei(t), 
we obtain that 



/(x + fct)(/)j(t)dt 



lim 

n — >oo 

= lim 



F„(x + fct)(?!)j(t)dt 



ai0i(x) / 0fei(t)0j(t)dt = 0. 
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If j = — fcl, then making a change of variables x + fct = s shows that 

/ /(x + fct)</>j(t)dt = / /(x + fct)(?!)i(-fct)dt = / /(s)0i(x-s)ds 
Jn Jn Jn 

which is a trigonometric polynomial in x in 7i„. Consequently, we conclude that 

PP'^ f is indeed a trigonometric polynomial in 7i„. 

Since J„.p(t) = J„_p(— t), it follows from (1) of Proposition and Minkowski's 

inequality that 



WFrf'fWp < 



K„Jt)Alfi-)dt 



< / if„.p(tK(/;|lt|l)pdt 
p Jn 



<c^,(/;i)p / KnAt){l+n\\t\\ydt<CLJrif;^)p, 
Jn 

where the last step follows from Lemma l473l □ 



For 1 < p < oo and r = 1, 2, . . ., define as the space of _ff-periodic functions 
whose r-th derivatives belong to L^. 

Corollary 4.5. For 1 < p < oo and r = 1, 2, . . ., i/ / G Wp then 

En{f)p<cn-'' II^'/IIp' n=l,2,.... 
|fe|=p 

4.4. Inverse Theorem. As in the classical approximation theory, the main task 
for proving an inverse theorem lies in the proof of a Bernstein inequality. For this, 
we introduce an operator that is of interest in its own right. Let 77 be a nonnegative 
C°° function on R such that 

ri{t) ==1, if < t < 1, and ri{t) = 0, if t < or t > 2. 

We then define an operator rjnf on Q by 

7?„/(x) / /(x- t)77„(t)dt, where 7?„(t) := ^ f - j i:)fe(t), 

where Dk is the Dirichlet kernel in 12.71 Evidently, rjnf G 'H2n and r]nf = / if 
/ € Tin- Such an operator has been used by many authors, starting from [8]. It 
is applicable for orthogonal expansion in many different settings; see, for example, 
[15j . A standard procedure of summation by parts leads to the fact that |l77n/|| 
is bounded. Consequently, using the fact that ijnf = f for / <E Hn, we have the 
following result. 

Proposition 4.6. For 1 < p < 00, if f ^ then 

WVnf - f\\p < cEM)p, n = l,2,.... 

This shows that for all practical purpose, rj^f is as good as the polynomial of 
best approximation. For our purpose, however, the more important fact is the 
following near exponential estimate of the kernel function ?7„(t). For a € Nq, write 
|a| = ai + . . . + ad- 

Lemma 4.7. For eaeh k ~ 1,2, . . there exists a constant Ck that depends on k, 
such that 

n\a\+2 _ 

a"ry„(t) < t e n. 
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Proof. The main tool of the proof is the Poisson summation formula ()2.1ip as used 
in the case of trigonometric series on the real line in [llj . Let us introduce a notation 
that \t\H = max{|ti|, \t2\, jial}- Since J(„ = {j e H„ : |ji| = n or \j2\ = n or jja] = n}, 
we can write r]„(t) as 



2n 



lJ|g 

n 



In particular, for a S Nq, we have 

a"^„(t) = (f)i"i^r;(^)r0j(t). 



jeH„ 



An immediate consequence of this expression is that 

(4.3) |9"ry„(t)| < cUU E ^ ^UU^^'^^^^ ■ 

jeH,. 

Define $„ such that $„(t) = 77 (l^jfi) t". Then $„(0 = /jja^ $„(t)e^« *dt. 
The definition of and Poisson summation formula shows that 



(4.4) 



77„(t) = 5»(j)'^j(t) = 273 ^ a>„(t + 3j). 



In order to estimate the right hand side we first derive an upper bound for <!>„. 
Since 77 is a function and is differentiable except when one of the variable 
is zero, taking derivatives in norm, we end up with 



;^)l''lt^$„(t)= / 9^$„(s)e^-Vs = / df' 



Nil I s" 



Each derivative of t] (^^^^ yields a n ^. For /3 G Nj] and k := |/3| > \a\, we have 



(^)t"] = E 



l7l<fe 

E 

l7l<|a| 



„-fc+l7l^(fc-|7l) 



„-fe+l7l^(fe-|7l) 



where we have used muti-index notations that for fc e N and a G Z'', a! = ai! . 
and (^) = fc!/(a!(fc-|a|)!). Hence, as ?7*^^^ is supported on [1, 2] for j > 1, we deduce 
that 



< en 



-I/3I + I"! 



^(fc-|7l)| 



E ii'? 

l7l<|a| 
fc 

J=fc-|a| 

Together with (|4.3p . we conclude then 

„|a|+2 



"<|t|-H<2n 



|$„(t)| < Ck 



(l + n||t||)'=' 



Cfc = c max |l?7*^^^|loo- 

h-\a\<3<k 
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As a result of the estimate, we conclude from (|4.4|) that 

|,„(t)| < c |<i>„(t + 3j)| < c. (l + ri|t'+3j||) 



k • 

~ \-L -|~ ~T~ '-'III/ 

Since ||t|| = inax{|ii|, Ita]} < 1 for t e rj, we have ||t + 3j|| > 3||j|| - 1 > 2|jj|| if 
j 7^ 0, and thus, 1 + n\\t + 3j|| > (1 + '^)||j||- Consequently 

and the proof is completed. □ 

Remark 4.1. The proof of the above estimate relies essentially on Poisson summa- 
tion formula, which is known to hold for all lattice AlA (see, for example, [71 llOj). 
Thus, there is a straightforward extension of the above result with an appropriate 
definition of partial sums. For relevant results on lattices, see [3]. 

As an application of this estimate, we can now establish the Bernstein inequality 
for hexagonal trigonometric polynomials. 

Theorem 4.8. If a ^ Nq, then for 1 < p < oo there is a constant Cp such that 

WSnWp < Cpnl"l|j^„||p, for all Sn G Hn- 
Proof. Recall that rjnf e 7i2n and rjnf ^ f for / e Tin- We have then 



^„(t) = (ry„5„)(t) = / Snis)rjn{t ~ s)ds 



For p — 1 and p = oo, we then apply the previous proposition with fc = 4 to obtain 

< cn^"\\Sn\\p f 77^7^4^*1^*2 < cnl"||5„||p, 

which establishes the stated inequality for p — I and p = oo. The case 1 < p < oo 
follows from the case of p = 1 and p = cx3 by interpolation. □ 

It is a a standard argument by now that the Bernstein inequality yields the 
inverse theorem. 

Theorem 4.9. There exists a constant Cr such that for each function f G C{n) 



UJ. 



0<n<h-^ 

5. Approximation on Triangle 

The hexagon is invariant under the reflection group A2 , generated by the reflec- 
tions in the edges of its three paris of parallel edges. In homogeneous coordinates, 
the three reflections (Ti, (72 and 173 are given by 

tax := -{ti,t3,t2), tcr2 := -(^2,^1,^3), taa -(^3, i2, ^i)- 

Indeed, for example, the reflection in the direction (\/3, 1) becomes reflection in 
the direction of a = (1, —2, 1) in Rj^, which is easy to see, using ti + t2 + ts — 0, 
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as given by t — 2j^^a — —(^2,^1, is) — io'2, The reflection group A2 is given by 

A.2 = {1, fJi, 0-2, CTs, 0-10-2, 0'2CTi}. 

Define operators V'^ and 'P~ acting on functions /(t) by 

(5.1) P±/(t) - i [/(t) + /(tai(72) + /(taacTi) ± /(tai) ± /(ta2) ± /(taa)] . 

They are projections from the class of _ff-periodic functions onto the class of invari- 
ant, respectively anti- invariant functions under A2- The action of these operators 
on elementary exponential functions was studied in and more recently studied 
in [131 m] and in [TD]. For 0k (t), we call the functions 

TCk(t) ■.= V+Mt), and TSk(t) := -V-^^) 

I 

a generalized cosine and a generalized sine, respectively. For invariant functions, 
we can translate the results over the regular hexagon to those over one of its six 
equilateral triangles. We choose the triangle as 

(5.2) A :={(ii, i2, h) : ii + i2 + is - 0, < ii, i2, -is < 1} 

={(ii,i2):ii,t2>0, h+t2<l}. 

It is known that the generalized cosines TCk are orthogonal with respect to the 
inner product 

(/,5)a := ^Jj{t)W)dt = 2 jj{h,t2)g{tiM)dtidt2, 

so that we can consider orthogonal expansions in terms of generalized cosine func- 
tions, 

/^^/kTCk, /k = (/,TCk)A, 

keA 

where A {k e H : fci > 0, A:2 > 0, A:3 < 0}. It is known that (/,.g>^ = (/,.9)a if 
fg is invariant. If / is A2 invariant and _ff-periodic, then it can be expanded into 
the generalized cosine series and it can be approximated from the space 

rc; = span{TCk : k e A, -/cg < n}. 

This is similar to the situation in classical Fourier series, in which even functions 
can be expanded in cosine series and approximated by polynomials of cosine. We 
state one theorem as an example. 

Theorem 5.1. /// G C(A) is the restriction of a A2 invariant Junction in C(il), 
then the (C,l) means of its generalized cosine series converge uniformly to f on A. 

We take this theorem as an example because the meaning of the (C, 1) means of 
the generalized cosine series should be clear and we do not need to introduce any new 
definition or notation. The requirement of / in the theorem may look redundant, 
but a moment of reflection shows that merely / S C'(A) is not enough. Indeed, 
in the classical Fourier analysis, an even function derived from even extension of 
a function / defined on [0, tt] (by f{—x) = f{x)) satisfies /(— tt) = /(tt), so that 
it is automatically a continuous 2tt periodic function if / is continuous. The H- 
periodicity, however, imposes a much stronger restriction on the function. Indeed, 
for a function / defined on A, we can then extend it to F defined on Q by A2 
symmetry. That is, we define 

F(t) = /(tcr), t e Act, a eA2. 
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It is evident that A = UaeA2'^'^- Ii^ order that is a continuous i?-periodic, we 
will need the restrictions of F on two opposite linear boundaries of Q are eqaul. 
Let dnA denote the part of boundary of A that is also a part of boundary of fl. 
Then 9siA := {((^1,^2,—!) : ^1,^2 > 0, ii + ^2 = 1}. Upon examining the explicit 
formula of <t G we see that F being continuous and iJ-periodic requires that 

F(t) =i^(tcr2), F{tai) ^ F{ta2ai), F{ta3) ^ F{taiCT2), t e OA. 
In terms of / this means 

(5.3) /(ti,i2,-l) = /(-t2,-ti,l), /(i2,-l,ti) 1,-^2), 

/(-l,ii,^2) = /(l,-i2,-ti), for ti+t2^l and h,t2 > 0. 

Hence, only functions that satisfy the restrictions (|5.3p can be extended to A2 
invariant functions on f2 that are also -ff-periodic and continuous on f2. As a result, 
we can replace the assumption on / in Theorem 1 5.1 1 bv *// € C(A) and f satisfies 
()5.3p . It does not seem to be easy to classify all functions that satisfy (|5.3|) , which 
is clearly satisfied if / has A2 symmetry. 
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